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POTENTIAL SOLAR BEAM IRRADIATION ON SLOPES 


Tables for 30° to 50° Latitude 


by 


Ernest C, Frank and Richard Lee 


Solar Energy Relations 


Solar irradiation, the energy source to the 
hydrologic cycle, photosynthetic processes, 
and other natural phenomena, must be evaluated 
continually in diverse land management and re- 
search ventures. The establishment and protec- 
tion of forests, and their yields of wood and 
water, are directly affected. Energy analyses 
for wild or agricultural lands invariably begin 
with data or assumptions relevant to the recur- 
ring input from the sun. In recent years this 


approach has emphasized the need for a quan- 
titative description of such concepts as 
"southern exposure," "solar climate," or "topo- 
graphic shading." Significant contributions 


have been those by Fons et al.(1960), Kaempfert 
and Morgen (1952), Loewe (1962), Milankovitch 
(1930), and Okanoue (1955). 


The slope and orientation of terrain ele- 
ments, especially important in forest hydrology, 
are important climatological factors as well. 
It is well known that the differential sunning 
of slopes alone can account for climatic ex- 
tremes measurable in terms of equivalent lati- 
tudinal differences. The experienced forester 
recognizes the extent of such differences in- 
tuitively, and local management practices may 
Ecrleet “his artiull tactics. Yet, scientific 
progress is facilitated to a greater extent by 
numerical qualification than by rules of thumb. 


Solar beam irradiation is only one of sever- 
al radiant energy components of importance in 
the radiation balance for natural areas, yet 
the differences between areas are often largely 
a function of this one parameter. Fairly uni- 
form climatic changes with latitude are 
obviously recognizable and related directly to 
sun angle. This relationship has been confirm- 
ed from long-term measurements of the various 


components 


of radiation exchange. The typical 


relation is: 


J+H+G-R-A (1) 


Wn 
Il 


where S is the net radiation balance, J the 
direct beam, H the diffuse sky, G the thermal 
atmospheric, R the reflected, and A the thermal 
terrestrial components. Data from two middle- 
latitude stations show that long-term averages 
of the additional terms H; G, R and A are 
largely self-canceling: 


Station (H + G)+(R + A) 
Davos, Switzerland 1959 0.961 
Hamburg, Germany 1957-62 1036 


Just as ordinary climatological data are not 
adequate to describe day-to-day events meteoro- 
logically, direct-beam irradiation alone cannot 
describe a short-term radiation budget. Con- 
versely, it is expected that mean values for 
this parameter will be as useful as mean tem- 
perature and precipitation values. 


Radiation detectors and integrators are 
being used more widely with increasing accept- 
ance of the usefulness of solar radiation 
values. This equipment is not used widely 
enough, however, to be a tool in the general 
classification of forest sites, empirical- 
theoretical radiation indexes have been evolved 
to meet such needs. Solar radiation values 
computed from the solar constant, an atmospher- 
ic attenuation constant (turbidity factor), and 
the geometry of a site with respect to the 
sun's rays, provide remarkably good long-term 
means. The kinds of parameters that change 
from day to day (clouds, smog, dust, and atmos- 
pheric moisture) make short-term forecasts much 
more uncertain. 


A purely theoretical parameter--potential 
solar beam trradtatton--is obtained by neglect- 
ing entirely the turbidity factor mentioned 
above (the influence of the atmosphere in de- 
pleting direct solar energy). The intensity on 
a surface is then proportional to the cosine of 
the angle of incidence. The incidence angle 
for a terrestrial surface depends in turn on 
the five independent variables: terrestrial 
latitude, hour angle (time of day), solar dec- 
lination (time of year), surface inclination 
or steepness, and surface orientation (azimuth 
of slope). 


solar beam irradiation 
appreciably from actual irra- 
diation, it is an appropriate, if nonrigorous, 
basis for comparing the energy flux among the 
facets of a landscape. In addition, it may be 
considered as a permanent site factor, the only 
variation being a perfectly cyclical one and 
calculable. The tables that form the bulk of 
this report are the result of such calculations. 
They provide a means of (1) quantitatively dis- 
tinguishing, in terms of potential insolation, 
between surfaces that differ as to latitude, 
slope inclination, and orientation, and (2) 
assessing the potential insolation received on 
selected dates by a specific surface, which 
facilitates seasonal comparisons. The tables 
give the potential solar beam irradiation in 
Langleys received during 1 solar day, and the 
sunrise-sunset times as hours before and after 
true solar noon for 24 dates during the year. 
Latitudes 30° and 50° North are covered by 2° 
increments. Slope inclinations from 0 to 100 
percent by 10-percent increments, and slope 
orientation to 16 points of the compass are 
covered. For inclined surfaces the latitude 
and change in longitude needed to locate the 
equivalent slope, or horizontal surface par- 
allel to the inclined surface, are tabulated. 
The total potential insolation received during 
1 year, and a radiation index (R.I.) based on 
the annual values, are given for all surfaces. 


Although potential 
usually differs 


Equations 


The mathematical theory underlying computa- 
tions of potential or extraterrestrial radia- 
tion was developed by Milankovitch (1930) in 
his classical 'Mathematische Klimalehre." In 
somewhat altered symbolism, the instantaneous 
irradiation Is of a horizontal surface is given 
as: 


Io 
Is = 2p (sin 6 + sin 6 
+ cos 6 * cos 5 * cos wt) (2) 
where: 
Io = solar constant 
e = radius vector, ratio of the 


earth-sun distance at a particular 
time to its mean 
= terrestrial latitude 
6 = solar declination 
angular velocity of the earth's 
rotation, 15 degrees per hour 


t = time in hours from noon 
The time of sunrise and sunset (ti> negative 
and t,, positive) is found by setting Is = 0 
and solving for t, or : 

Cos wt = -tan 6 « tan 6 (3) 

By integration of equation 2, the total 
radiant energy for the period of 1 solar day, 


Ig is found as 


Io 
ig, = 39 (Be °o gilim @ © sila © 


+ : 2 Cos 8) 0 Coe 9 20 silm we) (4) 

To account for the inclination and orienta- 
tion of sloping surfaces, Okanoue (1957) 
included the additional arguments, k = slope 
inclination, and h = slope azimuth. Slope in- 
clination, azimuth, and latitude are used to 
compute 6' and a as follows: 


‘ 


6' = sin ! (sin k * cos h * cos 6 
ar CoE 1 9 Sin ©) (5) 
a = tan ! [isin jn oO alin je) 
ad (OS Ik oO Com |) 
= COs} Je) O Spl Ig Oy Galion >) (6) 
These locate the equivalent slope (Lee 1962) 


which is the horizontal surface parallel to the 
inclined surface. The latitude of the equiva- 
lent slope is given by ®'. A negative 6" in- 
dicated an equivalent latitude south of the 
equator. The change in longitude (or hour 
angle) between the slope and its parallel hori- 


zontal surface is given by a. Okanoue's modi- 
fication of equation 2 with reference to the 
idea of an equivalent slope becomes: 
Is = io (Gin, @Y OO sin © 
e 
+ cos 6' * cos 6° cos wt") (7) 
where wt' = wt + a. Sunrise and sunset times 


for slopes are given by: 


cos wt' = -tan 6' + tan 6 (8) 
where it is understood that equation 3 contin- 
ues as the limiting case (i.e., day-length on a 
horizontal surface is the maximum for a given 
latitude; the inclined plane may have a shorter, 


but never a longer, solar day than the horizon- 


Table 1.--Input Data 


Solar 

Latitude Slope Slope inclination declination Dates, approx. ee Toler 
fe) : direction Deg. Min. Percent Deg. Min. 
30 00 0 00 0 ae 228) 30 6/22 0S 297, 1.93616 
32 00 N 5 43 10 ey 21 58 6/1 7/12 1.03090 1.94005 
34 00 NNE-NNW AH iS 20 ae) IG) 20 5/18 7/27 1.02728 1.94689 
36 00 NE=-NW 16 42 30 a ILS) Se) 5/3 8/10 1.02190 1.95714 
38 00 ENE-WNW Dall 48 40 + 10 35 4/19 8/25 LaOLa28 1.96990 
LO 00 E-W 26 34 50 + 5 38 4/4 9/9 1.00739 1.98553 
42 00 ESE-WSW 30 58 60 0 00 3/21 9/23 0.99960 2.00080 
44 00 SE-SW 35) 00 70 see) 38 S/T) 10/8 0.99154 2.01706 
46 00 SSE-SSW 38 40 80 - 10 55 2/20 10/22 0.98404 2.03243 
48 00 S 41 59 90 - 15 3)5) 2/7 11/5 0.97790 2.04520 
50 00 45 00 100 - 19 20 1/24 11/19 0.97285 2.05582 

=o 72IL 58 1/10 12/3 0.96938 2.06317 

=—23 30 12/22 0.96759 2.06699 


Values for the Radius Vector e interpolated after List (1958). 


tal facet at the same latitude). In certain 
instances double solar days may occur--both 
sunrise and sunset occur twice for the slope 
during one rotation period (24 hours). A 
common example is the north wall of a building 
in summer at a latitude (North) exceeding the 
declination in degrees. 


Over the time period of a day or less, the 
total potential slope irradiation, Iq, is de- 
rived from the relation: a 


Iq = = [ce2 (ey) alsa yO atau 
+ = co Cos GY oS CoS & 
(San wt,' - sin wy") | (9) 


1G Fe al oe , 
where ae is — in radians. 
W 


Explanation of Tables 


A Fortran program for the IBM 1620 computer, 
incorporating equations 3, 4, 5, 6, 8, and 9, 
was used to obtain the tabulated values. In- 
puts for this program were the Iote* values and 
the sine, cosine, and tangent (where necessary) 
of the latitudes, slopes, aspects, and declina- 
tions listed in table 1. The latitudes (30° 
to 50° North) include most of the conterminous 
United States, excluding only southern Texas 
and Louisiana, and most of Florida. While the 
potential insolation values listed apply equal- 
ly to slope aspects symmetric with respect to a 
north-south axis (NNE-NNW, E-W, etc.) the sun- 
rise and sunset times are for slopes with east- 
facing components only. Sunrise and sunset 


times on west-facing slopes are the reverse of 
those on east-facing slopes and, of course, 
their signs are reversed. In addition, the 
values of a for east-facing slopes are to the 
east, while for west-facing slopes the indi- 
cated longitude change is to the west. 


Interpretation of the tables for latitudes 
in the southern hemisphere” is straightforward 
if minor differences in season lengths are ig- 
nored. Slope aspects are reversed; there is a 
180° rotation of the north-south axis (N=S, 
NNE=SSW, etc.). The sign of the solar declina- 
tion is reversed, and the date changes by ap- 
proximately 6 months. Users of the tables can 
find corresponding dates in the southern hemi- 
sphere by inspection in table 1. There the 
corresponding date for the southern hemisphere 
is found by entering the table at the solar 
declination equivalent in value but opposite in 
sign to that for the date in the northern hemi- 
sphere. Sunrise and sunset times are given for 
southern hemisphere slopes with east-—facing 
components, just as for the northern hemisphere. 

True solar time was used exclusively in the 
computations. True solar time is reckoned as 
the number of hours’ before (-) or after (+) 
solar noon--the moment the geometric center of 
the sun crosses the meridian (longitude) of the 
observer. Likewise, for our purposes, the sun 
rises and sets when its geometric center 


3Little land area is included between 30° 
and 50° South Latttude--approximately the 
southern third of South America, southern Aus- 
tralia and New Zealand, and a very small part 
of South Africa. 


crosses the horizon. Sample calculations for 
converting the true solar times of sunrise and 
sunset used in the tables to local standard 


times are given in the appendix. 
The 13 solar declinations were chosen 
according to the approximate relation: 


) 


2BRo Suen aN (10) 
where N 
equinox. 
tain 
be 


is the number of days to the nearest 
Each quarter year is assumed to con- 

90 days, and solutions to equation 10 can 
arranged to yield declinations representa- 
tive of approximately equal numbers of days. 
Except for the solstice dates (June 22 and 
December 22) each declination is represented by 
two dates as shown in table l. 


The radius vector e varies seasonally, and 
is used as a correction for the solar constant 
to account for the varying earth-sun distance. 
Neglect of this correction introduces an error 
equal to the range of variation of solar energy 
outside the atmosphere, or about 6.7 percent. 
Lieywtacts that Ve ws net a simple Eunction oF 
declination 6 presents a minor complication. 
The procedure used was to take an arithmetic 
average of e for the 2 days represented by a 
single declination. The error involved in this 
approximation is less than one percent. 


Solutions to equations 5 and 6, representing 
the latitudes and changes in longitude of par- 
allel horizontal surfaces, are given in the 
tables for each slope-aspect combination. 
These will be interesting in a geographic 
sense, because they indicate to some degree the 


relative change in solar climate to be expected. 


Annual radiation totals based on a 365-day 
year, not 365.24 days, expressed in thousands 
Langleys (kg-cal/cm*) are derived by weighting 
each of the daily totals according to the num- 
ber of days for which each is representative. 
The R.I.'s for each slope are the ratios of 
annual radiation totals to the annual maximum 
potential solar beam irradiation. Since the 
mean day length at any latitude is 12 hours, 
the maximum potential insolation for the year 
is the solar constant times the duration of 
sunshine for the year: 


20) 365) *) 2) 760) = 525, 600 Langilieyse 


The indexes are directly comparable for all 
surfaces regardless of latitude. 


The R.I.'s, purely theoretical quantities, 
have application possibilities in forestry, 
ecology, and climatology as relative measures 
of the solar climate for sites. Land classifi- 
cation systems can well afford this elementary 
additional measure of the topoclimate. The 
indexes can also be used for more refined 
empirical-theoretical evaluations, as when ex- 
perience permits the assignment of a mean tur- 
bidity factor to a locality for a day, season 
or year. 
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Appendix: Sample Calculations 

Local standard time is based on a mean solar 
time for 15° longitudinal zones; Eastern Stand- 
ard Time, for example, is 75th Meridian time. 


The mean solar time for a locality is earlier 
(or later) then local standard time by 4 min- 
utes per degree of longitude measured east (or 


west) to the time zone meridian. Thus, mean 
solar time at Pittsburgh (80° West) is 20 min- 
utes earlier than Eastern Standard Time; at New 
Words Ghisy (A> Westie) te ag 4 imimnieas ilejecres 


solar time at a site differs from mean 
time through the approximate range, -14 
minutes to +16 minutes as shown in table 2. 
This difference is commonly referred to as the 
Equation of Time. From table 2, true solar 
noon on May 3rd occurs 3.1 minutes before mean 
solar noon. To complete the example for Pitts- 
burgh, the sun actually crosses the 80th merid- 
Lali Oil Way Shee aie LOO = Salk 16.9 minutes 
Buciceie “MOON, Ole Ze aloo@we IWS) joyoiily Weisieaiain 
Standard Time. And on February 7th, true solar 
NOOR dts News work City Occurs ate —4 4 1452 
10.2 minutes after noon, or at 12:10 EST. 


True 
solar 


Table 2.--Equation of Time, After List (1958) 
(Subtract s Equation, of Time from true solar 
time to obtain mean solar time). 


Jan 10 - 7.2 -0.12 
24 -12.0 -0.20 
Feb 7 -14.2 -0.24 
20 -13.9 -0.23 
Mar 7 -11.3 -0.19 
AL - 7.5 -0.13 
Apr 4 = 355) -0.05 
19 + 0.5 +0.01 
May 3 ap Si Ab +0.05 
18 + 3.7 +0.06 
June 1 + 2.4 +0.04 
22 - 1.6 -0.03 
Seley, 92 - 5.4 -0.09 
DY, - 6.4 -O.11 
Aug 10 - 5.4 -0.09 
25 - 2.3 -0.04 
Sept 9 ah Dod +0.04 
23 + 7.3 +05 12 
Oct 8 +12.1 +0.20 
22 = rls 3} +0.26 
Nov 5 +16.4 -Oe27/ 
19 +14.7 +O) 25 
Dec 3 ALO 5) =F Ole 
DD + 1.8 +0.03 


D 


In this way, any of the sunrise and sunset 
times in the tables can be adjusted to local 
standard time. For example: 


(a) The solar day for a 40 percent east- 
facing slope at Grand Canyon, Arizona (36° 
North, 112° West), on December 22nd begins at 
—4./7 hours and ends at +3.15 hours. The cor— 


responding Mountain Standard (105° Meridian) 


Times are found as follows: 


True solar times ANH Ines eral oP8\s Ib Invests 
Mean solar times SAO) IMesg! eyael  s-3jo 2 ines 
are earlier by 
OR03 phi. Geabilies2) 
Mountain Standard Aas omhics..e sand) | +359) sihiss). 
Times are later by 
I eS alOR ES aL. 
One PLS Ine, (O27) 
Mountain Standard U340) Acie — ermicl S}85)5) Daims 
Time 

(b) For a 40 percent west-facing slope at 
Grand Canyon on December 22nd, the solar day 
would begin at -3.15 hours and end at +4.77 
hours. The corresponding Mountain Standard 
Times are found as follows: 
True solar times 35 US lbeS 5 zbavel “SrA S 7/7) Inia) 
Mean solar times So INS} Ineo. Gievel sea 7/4) Inlies\e 
are earlier by 
OrO3. lwes  Cealoile 2)) 
Mountain Standard S2oYil Imeso. eiol GHG jwes, 
Times are later by 
LDF NOS) eS 
ore 7/5 Inee (Oss) 
Mountain Standard Qely Boils — Bal SSS). jo) iil 


Time 


(c) Sunrise and sunset March 2lst on a hori- 
VOM istic haces at Hort Collins. Colorado: (40. 
Noistin, WO = Wests) Cece ic) =O, 00 Eire 225,00) 
hours true solar time. The longitude coincides 
with the standard time zone meridian, so that 
the only correction required is the Equation of 
Time. 


True solar times =n 00Mhigs. and) —to.O00lhies. 
Mean solar times =H)58// Isese eiacl 44, 113} Ines. 
ane llater by 0.13 

Imes (esol 2)) 

Mountain Standard 6:08 a.m. and 6:08 p.m. 


Time 


THETA 


ALPHA 


SR SS 


LY 


YEAR 


List of symbols 


Equivalent latitude in degrees--the 
terrestrial latitude where a horizontal 
surface parallel or "equivalent" to the 
inclined surface is found. 


Change in longitude in degrees--the 
difference in longitude between the lo- 
cation of the inclined surface and its 
"equivalent" horizontal surface. Alpha 
is to the east for aspects with an east 
component, and to the west for aspects 
with a west component. : 


The elapsed time in hours from true 
solar moon at which the geometric 
center of the solar disc crosses the 
horizon. 


Total potential insolation in Langleys 
(gram calories/cm*) received by the 
surface during the solar day on the 
given date. 


Total annual potential insolation in 
thousands of Langleys (kilogram 
calories/cm*) received by the surface. 


Radiation index--the ratio of the total 
annual potential insolation to the max- 
imum potential insolation at the site. 
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